Entanglement and Quantum Phases in the Anisotropic Ferromagnetic Heisenberg 

Chain in the Presence of Domain Walls 
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We discuss entanglement in the spin- 1/2 anisotropic ferromagnetic Heisenberg chain in the pres- 
ence of a boundary magnetic field generating domain walls. By increasing the magnetic field, the 
model undergoes a first-order quantum phase transition from a ferromagnetic to a kink-type phase, 
which is associated to a jump in the content of entanglement available in the system. Above the 
critical point, pairwise entanglement is shown to be non-vanishing and independent of the bound- 
ary magnetic field for large chains. Based on this result, we provide an analytical expression for 
the entanglement between arbitrary spins. Moreover the effects of the quantum domains on the 
gapless region and for antiferromagnetic anisotropy are numerically analysed. Finally multiparticle 
entanglement properties are considered, from which we establish a characterization of the critical 
anisotropy separating the gapless regime from the kink- type phase. 
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I. INTRODUCTION 

In recent years, a great deal of effort has been de- 
voted to the understanding of entanglement in spin sys- 
tems. In fact, a number of schemes for quantum in- 
formation processing based on spin interactions have 
been proposed (e.g., Q, H H 0, 0|), motivating the the- 
oretical study of the amount of entanglement present 
in spin models USSUHmiHEElDl 
Furthermore, the interest in the content of entangle- 
ment available in statistical systems has also been in- 
creased due to the observed relationship between en- 
tanglement and the theory of qua ntum critical phenom- 
ena [H [13 H E HJ H |2l IH H. Indeed, it 
has been noted that quantum phase transitions present 
in condensed matter systems can be described from 
the point of view of entanglement. As examples, the 
ground state pairwise entanglement of XY and XXZ 
chains in a transverse magnetic field has been worked 
out [H EUll El UlEl, from which scaling behavior 
near the transition point is observed. For the transverse 
field Ising model, the quantum phase transition has also 
been characterized through the concept of localizable en- 
tanglement [23, which is the maximal amount of entan- 
glement that can be localized for two particles in the sys- 
tem by doing local measurements on the other particles. 
The typical scale for which the localizable entanglement 
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decays defines the entanglement l engt h, which diverges 
in the Ising quantum critical point |22J. Moreover, in the 
case of first-order quantum phase transitions, changes in 
the amount of pairwise entanglement at quantum critical 
points have been shown to be associated to macroscopic 
jumps in the magnetization for some frustrated spin mod- 
els \2',^ . In a general framework, the association between 
quantum phase transitions and bipartite entanglement 
has recently been analysed in Ref. |2^|- Hence, in ad- 
dition to its intrinsic relevance for quantum information 
applications, entanglement also plays an interesting role 
in the context of statistical mechanics. 

Keeping both motivations in mind, the aim of this work 
is to analyse the effect of a topologically non-trivial struc- 
ture on the entanglement of a spin system. More specifi- 
cally, we shall discuss the spin- 1/2 anisotropic ferromag- 
netic Heisenberg chain in the presence of domain walls. 
A related discussion about entanglement in spin models 
with open boundary conditions can be found in Ref. [2£| . 
The quantum domains can be generated by an antiparal- 
lcl magnetic field on the boundary spins of the chain. For 
a critical value of the magnetic field, the system under- 
goes a first-order quantum phase transition |2(3 |. which 
arises from the interplay among the exchange coupling 
J in the XY direction, the anisotropy A in the z-axis 
( A > J) , and the magnetic field h on the boundary. The 
quantum critical point separates a ferromagnetic phase 
from a kink-type phase, in which we show that entan- 
glement is present. It is interesting to observe that en- 
tanglement is known to be usually absent for ferromag- 
netic Heisenberg interactions regardless a magnetic field 
is present in a particular direction or the effect of temper- 
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ature is taken into account Even in the anisotropic 
chain, namely the XXZ model, entanglement is vanishing 
in the region A > J 0,0]. Actually it has recently been 
shown in Rcf. that entanglement can be nonvanishing 
for finite ferromagnetic Heisenberg chains at low temper- 
atures when a magnetic field is applied. However this 
entanglement is seen to disappear in the thermodynami- 
cal limit. The domain walls can thus be seen, due to the 
induction of the phase transition, as a general mechanism 
to promote the existence of entanglement in Heisenberg 
models with ferromagnetic couplings. Moreover, we will 
show that an analysis for large chains, in the kink-type 
phase, suggests that entanglement is independent of the 
boundary magnetic field. This result will allow for a gen- 
eral expression, in terms of the length of the chain, for 
the entanglement between arbitrary spins in the ground 
state, which is exactly known at the critical point. We 
shall also investigate the region A < J, corresponding to 
a gapless regime for — J < A < J and a gapful phase 
for A < — J. For these ranges of anisotropy, the distri- 
bution of nearest neighbor entanglement oscillates along 
the chain, exhibiting a clearly different behavior from the 
kink-type phase, where a monotonic increase is observed 
from the boundaries towards the center. Interestingly, 
a characterization of the critical anisotropy A = J can 
also be achieved from multiparticle entanglement prop- 
erties of the model. Indeed, we will show that A = ,/ 
corresponds to a minimum of the derivative of the global 
measure of entanglement introduced by Meyer and Wal- 
lach in Ref. jH- 

The paper is organized as follows. In the section II, 
we describe the quantum Heisenberg model for a ferro- 
magnetic domain wall, numerically discussing, for chains 
up to 24 sites, the general properties of the entanglement 
available in the system and its relation with the quantum 
phase transition. In Section III, we perform an analyti- 
cal computation of entanglement for general pairs of spins 
starting from the exact ground state wave functions at 
the critical point. Section IV is devoted to the numerical 
discussion of the case A < J. In Section V we analyse 
the entanglement properties of the model from the point 
of view of a multipartite measure. Finally, in section V, 
we summarize our results, presenting our conclusions. 



II. ANISOTROPIC HEISENBERG MODEL FOR 
FERROMAGNETIC DOMAIN WALLS 



h > represents the interactions of the spins with the 
boundary surfaces. The spin- 1/2 operators S", a — 
x,y,z, act on the site i (i = 1, ...,£) and are given by 
S a = cr Q /2, with a a denoting Pauli matrices. In order 
to focus the discussion on the values of the anisotropy and 
the magnetic field, let us assume, without loss of gener- 
ality, J = 1. As shown in Ref. [2g, the model presents a 
critical field 
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(2) 



that separates two quantum phases for a chain of arbi- 
trary length: a ferromagnetic ground state (h < h c ) and 
kink- type ground state (h > h c ). Remarkably, Eq. @ 
provides the exact location of the phase transition for 
chains of any size, the critical field remaining fixed as the 
number of sites is changed. Moreover, the ground state 
is exactly known at the quantum critical point, with de- 
generate wave functions given by j2(J 



I*o2)=^U^- JS3 >i, S2 ,..., S l>, (3) 
{«} 

where m = J^- Sf denotes the magnetization sector, Sj = 
±1/2, {|si, S2, sl)} is the basis where S z is diagonal, 
and M is the normalization factor, which is found to be 



-1/2 



(4) 



The sums over spin configurations {s} appearing in 
eqs. (0) and Q are to be taken in the sector of magne- 
tization m. The parameter q is defined by the following 
expressions 
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(q-q^ 1 ) and A = -(q + q- 1 ) 



(5) 



Before describing our results we present the definition 
of concurrence, which is the measure of entanglement 
used throughout this paper. As shown in Refs. [30L l3l|. 
there is a one-to-one correspondence between concur- 
rence and entanglement of formation [32j, the latter be- 
ing a strictly increasing function of the first. The con- 
currence Cp for a pair of qubits labelled as 1 and 2 is 
defined by 



The Hamiltonian for the one-dimensional Heisenberg 
model in the presence of a magnetic field generating do- 
main walls reads [26, 29] 



L-l 



L-l 



H = -J^2(SfSf +1 + SVSV +1 )-A^2SfS! +1 

i=l i=l 

-Hsi-si), 



(i) 



where the coupling J > and the anisotropy A > J 
are exchange parameters, and the effective magnetic field 



C12 = Max (Ai - A 2 - A 3 - A 4 , 0) , 



(6) 



where the Ai are the square roots, in decreasing order, of 
the eigenvalues of the operator 



(7) 



In Eq. J7J), pi2 denotes the density matrix, which can 
be either pure or mixed, for the pair of qubits 1 and 2, 
and p\ 2 hs complex conjugate in the the standard basis 
{| + +), | H — ), | — h), | }}. In a system with more 
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than two qubits, p\2 is obtained by tracing the density 
operator over the other qubits. The concurrence ranges 
from 0, implying an unentangled state, to 1, in which the 
two qubits are maximally entangled. 

As mentioned, the critical point given in Eq. (J2J is 
invariant under changes in the size of the chain. Actu- 
ally, as pointed out in Ref. (2(|, two sites are enough to 
identify the critical field. From the point of view of en- 
tanglement, the phase transition can indeed be observed 
by looking at jumps in the concurrence for arbitrary pairs 
of spins. For h < h c , the system is in the ferromagnetic 
phase, with vanishing concurrence for all the pairs. At 
the critical point, a wave function in the sector of mag- 
netization m — 0, which is entangled, becomes part of 
the ground state, but its mixture with the other sectors 
destroys entanglement. As h crosses the critical value, 
the ground state becomes non-degenerate, with compo- 
nents only in the sector m — 0, resulting in a jump 
in the concurrence for any pair of spins. This kind of 
non-analytical behavior in first-order transitions has been 
found in Ref. |23j |. where magnetization plateaus are as- 
sociated to jumps in the entanglement. In our case, nu- 
merical computation of concurrence up to 24 spins shows 
that the jumps are always observed in h — h c regardless 
the number L of sites, in agreement with Ref. f26|. 
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FIG. 1: Concurrence for the central pair of spins as a function 
of the number L of sites. The data have been obtained by 
taking the anisotropy A = 1.36 and the magnetic field h = 
0.5. As we can see, there is an exponential stabilization of the 
concurrence in the value C ~ 0.61, with the fitting function 
given by C = A + Ai exp(-A 2 L), where A = 0.612792, 
Ai = 0.359898 and A 2 = 0.669476. 

Concerning the distribution of entanglement along the 
chain, it is important to observe that the quantum do- 
mains break the ground state translation invariance, 
which implies that entanglement will depend not only 
on the distance of the spins, but also on their positions 
in the chain. The result is a concentration of pairwise 
entanglement in the central region, with the spins in the 
center forming the most entangled pairs. In fact, for a 
system with even number L of sites, it is possible to ob- 
tain a high degree of concurrence for the central pair, 



which evolves exponentially towards a stabilization value 
as we increase L. For instance, we show in Fig. 2] that 
this concurrence stabilizes in C « 0.61 for A = 1.36 and 
h = 0.5. This exponential behavior towards the stabi- 
lization is observed for all the pairs and the stabilization 
value is quickly obtained for large A and magnetic fields 
near h c . In fact, the stabilization value is strongly depen- 
dent on the anisotropy. In Fig. [2 we show how the con- 
currences for the pairs (L/2, L/2 + 1) and (L/2, L/2 + 2) 
change as a function of A in a chain with 22 sites. Note 
that, differently from the antiferromagnetic XXZ model 
without the quantum domains |20| . the maximum of C 
occurs for values of anisotropy greater than the isotropic 
value 1. 
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FIG. 2: Concurrence for the pairs (L/2, L/2 + 1) and 
(L/2, L/2 + 2) as a function of the anisotropy A. The data 
have been obtained by taking a chain with 22 sites and the 
magnetic field h = 5.0. The maximum value of C( £( / 2 ,l/2+i) 
occurs for A = 1.36 and of C(l/2,l/2+i) for A = 1.18. 

On the other hand, the study of the concurrence as a 
function of the magnetic field suggests that the stabiliza- 
tion value is independent of the field for large chains. For 
example, we show in Table[I]the concurrrence for the cen- 
tral pair as a function of L, for anisotropy A = 1.5 and 
several values of h. As it can be seen, there is clearly indi- 
cation to an approximately constant value C — 0.5974378 
for any value of h > h c assuming that we are not ex- 
tremely far from the critical field h c . As we will show 
in Section If f f I this numerical result leads to an analyti- 
cal derivation of a general expression for the concurrence 
between arbitrary spins since, as given by Eq. ©, the 
exact ground state of the model is known at the critical 
magnetic field for a general L. 



III. CONCURRENCE FROM THE GROUND 
STATE WAVE FUNCTION WITH VANISHING 
MAGNETIZATION 

In this section, we shall derive an analytical expression, 
holding in the kink-type phase, for the concurrences be- 
tween arbitrary spins along the chain. As will be shown, 
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L 


h = 0.6 


h = 5.0 


h = 20.0 


6 


0.599509933 


0.607452653 


0.607652241 


8 


0.597825866 


0.600881874 


0.601276968 


10 


0.597506824 


0.598224627 


0.598334941 


12 


0.597449660 


0.597593622 


0.597617282 


14 


0.597439774 


0.597466481 


0.597471030 


16 


0.597438106 


0.597442834 


0.597443658 


18 


0.597437830 


0.597438641 


0.597438785 


20 


0.597437785 


0.597437921 


0.597437945 


22 


0.597437778 


0.597437800 


0.597437804 


24 


0.597437776 


0.597437780 


0.597437781 



matrix, we have 



TABLE I: Concurrrence for the central pair C l / 2 ,l/2+i as a 
function of the chain size L for A = 1.5 and h = 0.6, 5.0, and 
20.0. Note that, as L is increased, Cl/2,l/2+i tends to the 
stabilization value 0.5974378 independently of h > h c . 



this formula comes from the analysis of the entanglement 
present in one of the degenerates ground states at h = h c , 
i.e. the wave function with magnetization m = 0. Thus, 
let us suppose that the system is prepared in the vanish- 
ing magnetization state |^o°i)' gi yen by Eq. ©. Then, 
the corresponding density operator is 



p T =N 2 9^ = lJSj+ ^ l=llS, |si.S2...SL}(Sis'2. 

{*}{*'} 



■s'l\ 



(8) 

where the sums over {s} and {s } are performed in the 
sector m = 0. The reduced density operator for a par- 
ticular pair of spins in the sites A and B can be written 
as 



P = Tr A,B( pT ^ 



(9) 



with Tr j B denoting trace over all degrees of freedom 
except the spins in the positions A and B. Therefore the 
matrix elements of p are 



= M 2 q A ( SA+s A) q B ( s B+s B ) x 

> jSj \ (10) 



P(SA,S B )\(S' A .S' B 



w 



where a a,, sb, s a , and s B represent the spin in the 
corresponding positions A and B. The sum over 
{s} stands for the set of (L — 2) spin configurations 
\siS2...S( A -i)S(A+i)—S(B-i)8(B+t)>»SL), whose magneti- 
zation is such that 



A-l 

E S J 

J=l 
A-l 

E s j 

j=i 



B-l 

- E -i 

3 =A+1 
B-l 

- E -i 



E -i 
E -i 



0, (11) 



0. (12) 



The conditions l|ll|l and (|12f) ensure that we are working 
in the sector of magnetization m — 0. Writing p as a 
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A - B n n o 
o o q 



-A-B 



with the functions fli (i = —1,0,1) given by 



(13) 



(14) 



where the sums over the (L — 2) spin configurations {si} 
are required to obey the restrictions 



'A-l 



B-l 



E=m5>+ £ Sj ' + -j +«=o. as) 

{s t } \j=l j=A+l j=B + l ) 

We are then able to compute the concurrence C, by 
means of Eq. (JBJ, for the spins in the sites A and B, 
which gives 



C = 27V 2 Q - \M fi-i 



(16) 



Hence we have obtained a simple and very general ex- 
pression to analyse concurrences in the model. In fact, 
as we discussed in Section [H] the ground state in the 
kink-type quantum phase is given by the wave function 
with m = and the corresponding concurrences are in- 
dependent of the magnetic field for large chains (see Ta- 
bic P). Therefore, concurrences obtained from the wave 
function with vanishing magnetization at h = h c , given 
by Eq. (|16fl . turn out to hold also in the kink-type phase, 
where h > h c . In other words, entanglement in a large 
chain, computed from the m — ground state wave func- 
tion, does not change as the magnetic field is increased. 
Eq. (|16f) displays a good property of being defined in 
terms of sums and not in terms of matrices, which are 
difficult to deal with when we have a large number of 
sites. It is worth mentioning that this expression holds 
as a good approximation even in the case of small num- 
ber L of sites assuming that the magnetic field is close 
to the critical value h c . Agreement with numerical simu- 
lations is observed, providing support to the exponential 
behavior, with respect to L, of the pair wise concurrences 
along the chain. 



IV. THE REGION A < 1 

As discussed above, the concurrence for any pair of 
spins, for A > 1, presents an exponential behavior to- 
wards a constant value in the kink-type phase. Moreover, 
the distribution of nearest neighbors concurrences along 
the chain exhibits a peak at the center, decaying as we get 
nearer the boundary. However, as A is changed to less 
than the critical value 1, entanglement properties become 
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rather different. As shown in Fig [31 the distribution of 
concurrences for nearest pairs are now characterized by 
an oscillating behavior, with a sequence of peaks along 
the chain. Observe that the oscillations become sharper 
as the anisotropy A is decreased and that the effect of 
the boundary field is a displacement of the peaks, with a 
damping of the concurrences. The peaks can physically 
be explained by the tendency to antiferromagnetic order- 
ing as A < 1, which is observed from the negative sign 
of the two-point correlation functions < S z (i) S z (i + 1) > 
for nearest neighbor spins i and i + 1. Thus, for h = 
the spins (1,2) tend to align antiferromagnetically, try- 
ing to form a bell state, and similarly the spins (3,4), 
(5,6), (L — 1,L) (for L even). The concurrences for 
these spins are greater than for the pairs (2, 3), (4, 5), 
(L — 2, L — 1) which, despite also exhibiting a tendency to 
antiferromagnetic alignment, are less correlated. When 
the domain walls are considered, the boundary spins are 
governed by the quantum domains, favoring the correla- 
tions of the pairs (2,3), (4,5), (L — 2,L — 1), which 
leads to a dislocation of the peaks, in addition to a reduc- 
tion of their sharpness. In the next section, we will show 
that the critical anisotropy A = 1 can also be very well 
characterized by looking at multiparticle entanglement 
properties of the model. 
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FIG. 3: Distribution of entanglement for A < 1 in a chain 
of L = 22 sites. C(i) represents the concurrence between the 
sites i and i + 1, where i = 1, L — 1. The long-dashed lines 
with circles are associated to ft = and the solid lines with 
squares to h — 3.0. In (a) we show the results for the gapless 
region for A = 0.5 and in (b) concurrence is plotted in the 
gapful phase for A = —1.5. 



V. MULTIPARTITE ENTANGLEMENT 

Let us discuss now some multiparticle entanglement 
properties of the model, which will be seen to provide 



support to the results presented before. It is impor- 
tant to emphasize that, for multiparticle systems, we 
do not have a unique measure of entanglement. Then 
it turns out that a number of measures have been pro- 
posed |33|, yj, |35l Ha, |33, |38|, |33 and different mea- 
sures can be related to distinct physical aspects of mul- 
tipartite entanglement. We shall consider here the so- 
called global entanglement, which is a scalable measure 
for pure states introduced by Meyer and Wallach in 
Ref. [2^|. For a general quantum state \tp) in the Hilbcrt 
space (C 2 )® L , the global entanglement is 

A L 

Q(i^)) = z E^(i u0) )'i w0) ))' 
3=1 

where the vectors and \v^}, which are elements of 

the Hilbert space (C 2 )® i_1 , are defined by 

2 L-1 

|«W))=t i (-)|V»)= Y.v&lx), 

05=1 
2 L-1 

\ v W)=t j (+M='£v(j»\y). (18) 
y=i 

with tj(b), for b 6 {+,—}, denoting the following opera- 
tion over basis vectors |6i&2---&l) of (C 2 )® L : 

t j {b)\b 1 b 2 ...b L ) = 8 bhj \b 1 ...b j ...b L ), (19) 

In Eq. I|19[) . the symbol * means absence. The function 
D(\uW),\vW)) in Eq. (Q2J is the norm-squared of the 
wedge product of \u^) and \v^'} 

D = E |u£*>t#) _ u««W| 2 . (20) 

x<y 

It has recently been shown in Ref. 01 that the global 
entanglement can be expressed in a simple way in terms 
of the one-qubit reduced density operator of the system 

Q(m=2^1-±pTr(p 2 J ) ) j , (21) 

where pj is the density matrix for the spin j after tracing 
out the rest. From Eq. 1)21(1 we can numerically compute 
the global entanglement in our model. Since the measure 
is defined only for pure states, we shall always be referring 
to the kink- type phase, where h > h c . As a first property, 
and in agreement on the results for the concurrence, Q 
is seen to be independent of the magnetic field if we are 
not too far from the critical point. This motivates the 
analytical computation, as performed in the case of the 
concurrence in Section IIIII of the global entanglement 
from the sector m = of the exact ground state Eq. JjJJ • 
The one-qubit reduced density matrix for a site A in this 
case reads 
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where N is the normalization factor given by Eq. J2J and 
the functions fli (i = —1/2, 1/2) are defined by 

fi. = J2 gKD^i 1 ^ + X£=a + i i'O , (23) 
Kl 

with the sums over the (L — 1) spin configurations {si} 
obeying the restrictions 

E=mX> + E s H +i=0 - ( 24 ) 

Thus the global entanglement Q can analytically be ex- 
pressed by 

The expression above can be used to compute, in a very 
good approximation, the global entanglement for large 
chains in the kink-type phase assuming that we are not 
extremely far from the critical point, since a very high 
antiparallel boundary magnetic field can considerably in- 
crease Q. 
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Fig-El This characterization of the transition is rather at- 
tractive since, based on a multipartite measure of entan- 
glement, it reflects a qualitative change in the variations 
of the ground state wave function of the system. Related 
characterizations of quantum phase transitions in terms 
of the global entanglement can be found in the Ref. |4l| 
for the transverse field Ising model and in Ref. ^2 f° r 
the single-mode Dicke Hamiltonian. 

VI. CONCLUSION 

In this paper, we have presented a discussion of how 
a nontrivial geometric interface, the ferromagnetic do- 
main wall, can modify the pattern of entanglement in 
the spin- 1/2 anisotropic Heisenberg chain. As shown, the 
quantum domains can be seen as a mechanism to gener- 
ate pairwise entanglement in a model with Heisenberg 
ferromagnetic couplings, which is usually uncntanglcd. 
Moreover we have seen that the amount of entanglement 
generated can be relatively large, specially for the central 
pair of spins. These effects are consequence of the first- 
order quantum phase transition induced by the presence 
of the domain walls. As a further nice result, we have 
been able to compute an analytical general expression for 
the concurrence of large chains thanks to the numerical 
observation that, in that case, entanglement is indepen- 
dent of the boundary magnetic field. Finally, we have 
detected different patterns for the entanglement in the 
kink-type and A < 1 regimes. From the point of view 
of pairwise entanglement, these regions are identified by 
a different distribution of entanglement along the chain. 
Moreover these distinct phases have also been shown to 
be characterized by looking at the problem from the point 
of view of multiparticle entanglement. By considering the 
derivative of the global entanglement with respect to the 
anisotropy we showed that a minimum occurs exactly at 
the critical point A = 1. Applications of this kind of 
phase characterization in other models and the possible 
role of other multipartite entanglement measures are left 
as open problems for future investigation. 



FIG. 4: Derivative of the global entanglement Q with respect 
to the anisotropy A. The magnetic field is set to h = 1.0. 
Note the minimum in the critical value A = 1. 

A further nice porperty of the global entanglement is 
that it presents an inflection point which characterizes 
the transition from the kink-type phase to the regime 
A < 1. In fact computing, from Eq. I|21|). the derivative 
of Q with respect to the anisotropy, at a fixed number of 
sites, we observe the existence of a minimum exactly in 
A = 1, which gets more pronounced as we increase the 
number of sites. This behavior is explicitly displayed in 
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